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ABSTRACT 

We present Very Large Array observations of twelve late M and L dwarfs in the Solar neighborhood. The observed 
sources were chosen to cover a wide range of physical characteristics — spectral type, rotation, age, binarity, and 
X-ray and Ha activity — to determine the role of these properties in the production of radio emission, and hence 
magnetic fields. Three of the twelve sources, TVLM 513-46546, 2MASS J0036159+182110, and BRI 0021-0214, 
were observed to flare and also exhibit persistent emission, indicating that magnetic activity is not quenched at the 
bottom of the main sequence. The radio emission extends to spectral type L3.5, and there is no apparent decrease in 
the ratio of flaring luminosities to bolometric luminosities between M8-L3.5. Moreover, contrary to the significant 
drop in persistent Ha activity beyond spectral type M7, the persistent radio activity appears to steadily increase 
between M3-L3.5. Similarly, the radio emission from BRI 0021-0214 violates the phenomenological relations 
between the radio and X-ray luminosities of coronally active stars, hinting that radio and X-ray activity are also 
uncorrected at the bottom of the main sequence; an even stronger violation was found for the brown dwarf LP944- 
20. The radio active sources that have measured rotational velocities are rapid rotators, vsini > 30 km sec"', while 
the upper limits on radio activity in slowly -rotating late M dwarfs (vsini < 10 km sec"') from this survey and from 
the literature are lower than these detections. These observations provide tantalizing evidence that rapidly-rotating 
late M and L dwarfs are more likely to be radio active. This possible correlation is puzzling given that the observed 
radio emission requires sustained magnetic fields of ^ 10- 10^ G and densities of ^ 10'^ cm"^, indicating that the 
active sources should have slowed down considerably due to magnetic braking. 

Subject headings: stars: low mass,brown dwarf-stars: activity-stars: radio emission-stars: magnetic 
fields-radiation mechanisms: nonthermal 



1. INTRODUCTION 

Activity in dwarf stars, as measured through the ratios of 
Ha and X-ray luminosities to bolometric luminosities, drops 
significantly in spectral types later than M7, and is possibly 
dominated by flai-es (Liebert et al. 1999; Reid et al. 1999; 
Gizis et al. 2000; Kirkpatiick et al. 2000). In addition, the 
dependence of activity on rotation appears to break down in 
sources with spectral type later than early-mid M, so that even 
rapidly rotating sources (vsin/ > 20 km sec"') show a decrease 
in emission (Basri & Marcy 1995; Tinney & Reid 1998) . On 
the other hand, the transition from a radiative-convective (ail) 
to a turbulent (a^) dynamo (Durney, De Young & Roxburgh 
1993) around the convective mass limit of 0.3 Mq (spectral 
type ^ M3) is not accompanied by a drop in activity (Hawley, 
Gizis & Reid 1996), indicating that in early M dwarfs the turbu- 
lent dynamo can sustain active coronae. It has been suggested 
that the observed shift to weak flaring activity in late M and L 
dwarfs is possibly due to a change in the nature of the dynamo, 
or decreased photospheric ionization levels. If radio activity 
follows the same pattern, then late M and L dwarfs would ex- 
hibit at most only weak flaring, and possibly weaker quiescent 
radio emission. 

In addition, Guedel & Benz (1993), and Benz & Guedel 
(1994) have shown that for a broad range of coronally ac- 
tive stars (up to spectral type M7) there is a correlation be- 
tween the radio and X-ray luminosities. Simply, the quies- 
cent luminosities are r elated by Lrad.q « Lx,q/10'^-^ Hz"', while 
the relation for flaring luminosities is slightly non-linear with 



Lrad,! ~ ^x.f/lO'^^ Hz ' for radio flares with L^-^d ~ lO'^ erg 
sec"' Hz"'. These relations are explained in terms of a heat- 
ing process in which there is a causal relation between the 
radio-emitting electrons and the X-ray emitting thermal plasma. 
The extension of the Guedel-Benz relations to spectral types 
later than M7 based on X-ray observations (Neuhauser et al. 
1999; Rutledge et al. 2000) predicts radio emission well below 
the detection limit of the Very Large Array. Past observations 
seemed to confirm this prediction (Krishnamurthi, Leto & Lin- 
sky 1999). 

The recent detections of an X-ray flare (Rutledge et al. 2000) 
and persistent and flaring radio emission (Berger et al. 2001; 
hereafter BOl) from the brown dwarf LP 944-20 (spectral type 
M9) pose a serious challenge to these ideas, since they demon- 
strate that late dwarfs can exhibit strong activity closely related 
to magnetic fields. The quiescent radio emission in particular 
indicates that the magnetic activity is persistent, contrary to Ha 
and X-ray observations of this source (Rutledge et al. 2000). 
Moreover, the observed emission violates both Guedel-Benz 
relations by approximately four orders of magnitude despite the 
fact that LP 944-20 is only an M9 source, hinting at a possible 
breakdown in the relation between radio and X-ray activity be- 
tween spectral types M7 and M9. Finally, the rapid rotation of 
LP 944-20, vsin/ « 30 km sec"', is puzzling given that the radio 
emission requires the presence of magnetic fields and ionized 
material, which should lead to magnetic braking. 

These observations raise several crucial questions: Is LP 944- 
20 a unique object, or do other late-type stars and brown dwarfs 
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show similar levels of radio activity? What is the underlying 
mechanism that gives rise to the emission, and is it similar to 
the mechanism in early M dwarfs? Is the same process respon- 
sible for the observed Ha and X-ray flares, and in some cases 
persistent emission? Is the activity correlated with particular 
physical properties? 

To address these questions, we have undertaken a survey of 
twelve late-type dwarfs (ranging from M5 through T6) that span 
a wide range of properties. The purpose of this survey is to 
search for flaring and/or persistent radio emission, and to in- 
vestigate whether this emission, and presumably magnetic field 
production, are correlated with particular physical properties of 
the sources. 

We summarize the properties of the survey sources in §2 
The radio observations and data reduction are detailed in §3 
In §Q we show that three of the twelve sources exhibit radio 
emission. We analyze the radio emission, show that at least in 
one case it violates the Guedel-Benz relations, and estimate the 
magnetic field strengths and coronal densities in §|[ We also 
discuss the implications for coronal heating mechanisms (§ 5.3 ). 
Finally, in §|| we compare the results of the survey to previous 
radio surveys of M and L dwarfs, as well as LP944-20, and we 
draw preliminary conclusions about possible correlations be- 
tween radio activity and physical parameters. 

2. TARGET SELECTION 

We observed twelve late M and L dwarfs, two of which are 
confirmed brown dwarfs, ranging in distance from « 8 - 160 
pc, and spanning a wide range of ages, rotation, activity, and 
spectral types. Below and in Table |l] we summarize the main 
properties of each source in the survey. 

CRBR 15: This source is a brown dwarf candidate in the p Oph 
molecular cloud (Wilking, Greene & Meyer 1999). The mass 
and age of CRBR 15, M < 0. 1 M© and f < 3 Myr, are estimated 
by a comparison to evolutionary models and isochrones (Bur- 
rows et al. 1997; Wilking, Greene & Meyer 1999). Wilking et 
al. (1999) claim that while the age of this source is secure rel- 
ative to other sources in the p Oph cloud, the absolute scale is 
model-dependent. Still, this source is the youngest in our sam- 
ple. The temperature and bolometric luminosity of CRBR 15 
are estimated at 2900-3000 K and log(Lboi/i0) « -1.0 to 
-1.3 (Luhman & Rieke 1999; Wilking, Greene & Meyer 1999). 
In a recent observation of the p Oph molecular cloud with 
Chandra (Imanishi, Tsujimoto & Koyama 2001) an upper limit 
of Lx < 1-1 X 10^^ erg sec"' was found for this source, with 
log(Lx/iboi) < -4.2, well below the X-ray saturation limit, 
log(Lx/iboi) ~ -3. In addition, CRBR 15 exhibits excess emis- 
sion at 2.2 pm (Wilking, Greene & Meyer 1999), and has a 
flux density of 40 mJy at 1.3 mm (Motte, Andre & Neri 1998). 
These observations indicate that it is associated with a circum- 
stellar envelope or an accretion disk. Thus, radio emission from 
this source can be due to accretion or deuterium burning. 
LHS 2243: This source is a slowly rotating, vsin/ < 5 km sec"', 
M8 dwarf. It shows evidence for weak persistent Ha emis- 
sion, as well as Ha flares, with an increase of ^ 30 in the line 
equivalent width (Martin, Rebolo & Magazzu 1994; Gizis et al. 
2000). The temperature of LHS 2243 is estimated to be - 2900 
K (Martin, Rebolo & Magazzu 1994) based on the calibration 
of Kirkpatrick et al. (1993). This calibration provides the high- 
est temperatures among the published temperature conversions 
(Martin, Rebolo & Magazzu 1994). Bessell (1991) suggested 



that this source is possibly as young as a few x 10^ yr, based on 
its smaller proper motion compared to bluer stars. For an age 
of ~ lO** yr, the inferred mass, based on an upper limit on Li 
abundance and several evolutionary models is ^ 0.06-0.07 M© 
(Martin, Rebolo & Magazzu 1994). If the source is younger 
than 10'' yr, the mass would be slightly higher. 
TVLM 513-46546: This is a nearby M8.5 dwarf (Tinney 1993; 
Tinney et al. 1995; Reid et al. 2001). TVLM513-46546 
(hereafter TVLM 5 13) is the fastest rotator in our sample with 
vsim w 60 km sec"', but it has only weak Ha emission with 
equivalent width of 2.5A (Martin, Rebolo & Magazzu 1994). 
Based on 1-2.5 pm spectroscopy, and synthetic spectra fitting, 
Leggett et al. (2001) find log(Lboi/Lo) « -3.65 and T w 2200 
K for this source. We note however, that their temperature scale 
provides consistently lower temperatures relative to other esti- 
mates (see notes on LHS 2065 below). 

LHS 2065: This source has spectral type M9 and relatively 
slow rotational velocity, vsin/ « 9 km sec"' (Tinney & Reid 
1998; Martin et al. 1999; Reid et al. 2001). It is active in 
Ha, exhibiting strong flares with a recurrence of < 0.03 hr"', 
and possibly more frequent weak flares at a rate ^ 0.5 hr"' 
(Martin & Ardila 2001). There are several estimates of the 
effective temperature of this source in the literature, ranging 
from 2100-2600 K. Leggett et al. (2001) provide an esti- 
mate of 2100 K (see notes on TVLM513 above). Another es- 
timate, based on near infrared photometry, provides a temper- 
ature of 2300 K (Leggett, Aflard & Hauschildt 1998; Martin 
et al. 1999), while Martin, Rebolo & Magazzu (1994) pro- 
vide a value of 2630 K (see notes on LHS 2243 above). The 
luminosity of LHS 2065 is estimated to be log(L/Lboi) ~ -3.5 
(Leggett et al. 2001), and its mass is estimated at ^ 0.06-0.1 
Mq (Martin, Rebolo & Magazzu 1994) which places it close to 
the brown dwarf cutoff. 

BRI 0021-0214: This nearby, M9.5 dwarf is a rapid rota- 
tor, vsin/ K, 40 km sec"' (Basri & Marcy 1995; Tinney & 
Reid 1998), but it exhibits no persistent Ha emission (Tin- 
ney & Reid 1998). Reid et al. (1999) discovered an Ha 
flare from this source, with a peak emission three times lower 
than the mean quiescent emission of early M dwarfs. In addi- 
tion, BRI 0021-0214 (hereafter BRI 0021) has an upper limit 
on X-ray emission of \og{Lx/Lbo\) < -4.7, well below the X- 
ray saturation limit of log(Lx/i'boi) ~ -3 for X-ray active stars 
(Neuhiiuser et al. 1999). Leggett et al. (2001) estimate the 
temperature of this source at 2100 K, and its luminosity at 
log(Lboi/i'o) ~ -3.4. Finally, Martin, Basri & Zapatero Oso- 
rio (1999) find an upper limit on Li EW of < 0.08 A. The rapid 
rotation, flaring emission, and lack of persistent Ha and X-ray 
emission indicate that this source is an analog of LP 944-20. 
PC 0025H-0447: This source is possibly a young (< 1 Gyr) 
brown dwarf (M< 0.06 Mq; Martin, Basri, & Zapatero Os- 
orio 1999), based on the detection of Li with EW= 1.0 ±0.3 
A. PC 0025H-0447 (hereafter PC 0025) exhibits very strong and 
persistent Ha emission (Graham et al. 1992; Mould et al. 
1994). Burgasser et al. (2000), on the other hand, suggest 
that PC 0025 is an interacting binary in order to explain the 
Li and Ha observations. The rotational velocity of PC 0025 
is surprisingly low, vsin/ « 13 ± 3 km sec"', given the strong 
level of Ha activity. In addition, PC 0025 has an upper limit 
on X-ray emission of logLx < 27.2 erg sec"' (Neuhauser et 
al. 1999). The bolometric luminosity of this source is esti- 
mated at log(Lboi/i0) ~ -3.7 (Schneider et al. 1991), giving 
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log(Lx/i'boi) < -2.7. The slow rotation and strong Ha activ- 
ity place PC 0025 on the opposite end of the activity-rotation 
spectrum from BRI 0021 and LP 944-20. 
2MASP J0345432H-254023: This source is an isolated LO 
dwarf, but it is still not clear whether it is a low mass star or 
a brown dwarf younger than ^ 1 Gyr (Kirkpatrick, Beichman 
& Skrutskie 1997). Leggett et al. (2001) estimate the luminos- 
ity of 2MASP J0345432H-254023 (hereafter 2MASP0345H-25) 
at log(Lboi/i'o) ~ -3.6, with an effective temperature of 2000 
K (see notes on TVLM513 above). 2MASP 0345H-25 has an 
upper limit on Ha emission of log(LHa/i'boi) <, -5.6. 
2MASSW J0036159H-182110: This source is a recently- 
discovered L3.5 candidate brown dwarf (Kirkpatrick et al. 
2000; Gizis et al. 2000; Reid et al. 2000). 2MASSW 
J0036159H-182110 (hereafter 2MASS 0036H-18) shows no ev- 
idence of Li or Ha, with upper limits of EW< 0.1 A (Reid 
et al. 2000; Kirkpatrick et al. 2000). The temperature of 
2MASS 0036H-18 is estimated to be 1800 K based on the spectra 
and models of Leggett et al. (2001). 

GD 165B: This source is a rapidly-rotating, vsini w 37 km 
sec"', L4 dwarf orbiting a white dwarf primary (Becklin & 
Zuckerman 1988). Kirkpatrick et al. (1999) used model atmo- 
spheres, including the effects of condensation and dust opac- 
ities, to estimate the temperature of this source, T w 1900 K, 
and in addition used the cooling age of the white dwarf pri- 
mary to estimate the age of GD 165B at 1.2-5.5 Gyr Based 
on these two quantities they conclude that this source is likely 
a brown dwarf. The luminosity of GD 165B is estimated at 
log(Lboi/i0) « -4.1 (Leggett et al. 2001). Along with DENIS- 
P J1228.2-1547 (see below), and possibly PC 0025, observa- 
tions of GD 165B can be used to test the effect of binarity on 
radio emission from late dwarfs, and the 4" separation between 
the binary members allows us to pinpoint the source of any ra- 
dio emission from the system. 

2MASSW J1507476-162738: This source is a nearby L5 
dwarf (Reid et al. 2000,2001). 2MASSW J1507476-162738 
(hereafter 2MASS 1507-16) has only upper limits on Li and 
Ha of < 0.1 A and < 0.5 A, respectively (Kirkpatrick et al. 
2000; Reid et al. 2000). 

DENIS-P J1228.2-1547: This source is an L5 brown dwarf 
exhibiting relatively slow rotation, vsin/ « 11 km sec"', and 
no detectable Ha emission, log(L//c(/Lboi) < -5.3 (Martin et 
al. 1997; Delfosse et al. 1997; Tinney, Delfosse & Forveille 
1997). The bolometric luminosity of DENIS-P J1228.2-1547 
(hereafter DENIS 1228) is estimated to be log(Lboi/i0) ~ -4 
to -4.3 (Leggett et al. 2001; Martin et al. 1997; Delfosse et 
al. 1997; Tinney, Delfosse & Forveille 1997), and its effective 
temperature is 1800 K (Leggett et al. 2001). In addition, DE- 
NIS 1228 has an upper limit on X-ray emission of logLx < 27.3 
erg sec"', and log(Lx/Lboi) < -2.3. Based on the detection of 
Li and the estimated effective temperature, the age and mass of 
DENIS 1228 are estimated to be M< 0.065 Mq and f < 1 Gyr 
(Tinney, Delfosse & Forveille 1997). Finally, Martin, Brandner 
& Basri (1999) have shown that DENIS 1228 is in fact a brown 
dwarf binary with a projected separation of approximately 5 
AU. While our observations cannot resolve the system, we can 
study the combined emission, and effect of binarity on the pro- 
duction of this emission. 

SDSS J134646.45-003150.4: This is the only methane (spec- 

'The VLA is operated by the National Radio Astronomy Observatory, a fi 
Associated Universities, Inc. 



tral class T6; Burgasser et al. 2001; Geballe et al. 2001) 
brown dwarf in the sample. It has a slightly higher effec- 
tive temperature than that of two other known methane brown 
dwai-fs, GL229B and SDSS J162414.37H-002915.6 (Tsvetanov 
et al. 2000), and it is inactive in Ha, with log(LHQ/iboi) <; -5.3 
(Burgasser et al. 2000). The bolometric luminosity of SDSS 
J 1 34646.45-003 150.4 (hereafter SDSS 1 346) is log(Lboi /Lq) « 
-5.3 (Burgasser et al. 2001; Tsvetanov et al. 2000). At a dis- 
tance of only - 11 pc (Tsvetanov et al. 2000) SDSS 1346 is 
an excellent candidate for the detection of radio emission from 
this class of objects. 

3. RADIO OBSERVATIONS 

Very Large Array (VLA') observations were conducted at 
8.46 GHz in the standard continuum mode with 2 x 50 MHz 
contiguous bands. A log of all observations is given in Table |[ 
Following the detection of emission from BRI 0021 at 8.46 
GHz, we also observed this source at 4.86 GHz. In all obser- 
vations we used the extra-galactic sources 3C48 (J0137H-331) 
and 3C286 (J1331H-305) for flux calibration, while the phase 
was monitored using calibrator sources within ^ 5° of the sur- 
vey sources. 

The first set of observations was undertaken in the A and 
A^B configurations. Unfortunately, half the data in these ob- 
servations had to be discarded due to a malfunction of the lower 
Intermediate Frequency channel, and the resulting rms noise 
levels are therefore higher than the theoretical values by w \/2. 
We repeated a subset of these observations over a period of ap- 
proximately three weeks in the CnB configuration. The final 
set of observations, with five new targets was undertaken in the 
DnC configuration. 

The data were reduced and analyzed using the Astronomical 
Image Processing System (AIPS, release 31DEC1999; Foma- 
lont 1981). Both the visibility data and maps were inspected 
for data quality and noisy points were removed. To search for 
flares, we constructed lightcurves using the following method. 
We removed all the bright background sources in each image 
using the AIPS/IMAGR routine to CLEAN the region around 
each source, and the AIPS/UVSUB routine to subtract the re- 
sulting source models from the visibility data. We then plotted 
the real part of the complex visibilities at the position of each 
target source as a function of time using the AIPS/UVPLT rou- 
tine. The background source subtraction is necessary since the 
side-lobes of these sources, and the change in the synthesized 
beam shape during the observation result in flux variations over 
the map, which may contaminate any real variability or gener- 
ate false variability. 

The uncertainties in the resulting lightcurves (Figures |l]-^ 
are dominated by system thermal noise. The system tempera- 
ture includes contribution from background sources, the sky, 
and hardware (e.g. receivers). At the observing frequency 
used in this survey, 8.46 GHz, the dominating term is the re- 
ceiver noise. As a result, there are no significant systematic ef- 
fects since in all cases the sources were located within a few 
arcseconds of the array pointing center, and the background 
source subtraction successfully removed possible contaminat- 
ing sources. In fact, the plotted uncertainties, as well as the rms 
noise from the maps, are within a few percent of the theoretical 

of the National Science Foundation operated under cooperative agreement by 
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noise estimates for the VLA at 8.46 GHz. 

To further test the veracity of the resuhing Hghtcurves and 
uncertainty estimates we constructed equivalent Hghtcurves for 
random positions on each map. We find that on average the 
fluctuations in these Hghtcurves do not exceed 2a, and the noise 
estimates are similar. 

4. RESULTS 

We detected flares and persistent emission from TVLM513, 
2MASS 0036+18, and, with lower significance, from BRI 0021. 
Figures |l|-Q show the Hghtcurves for observations of these 
sources. The peak flux densities range from w 360 /iJy in 
BRI 0021 to w 980 /iJy in TVLM513, and the flare timescales 
(FWHM) range from w 6 min in BRI 0021 to w 20 min in 
2MASS 0036+18. The resulting duty cycles, defined as the ra- 
tio of the time the sources produce flares to the total observing 
time, are « 2 - 10%. In the case of BRI 0021, the value of 2% 
is consistent with a duty cycle of < 7% for Ha flares (Reid et 
al. 1999). The duty cycles of all three sources are consistent 
with an average Ha flare duty cycle of ~ 7% for late M dwarfs 
(Gizis et al. 2000). 

The fraction of circular polarization during the flares is sig- 
nificant in TVLM513, /ci,c ~ 66 ± 4%, and 2MASS 0036+18, 
/circ ~ -62 ±5%, but is somewhat lower in BRI 0021, /cue ~ 
30 ± 20% (see Figures [l| and |[ and Table |]). The negative value 
for 2MASS 0036+18 indicates that the emission is left-handed 
circularly polarized. 

In addition, we clearly detect persistent emission from 
TVLM513 and 2MASS 0036+18, with flux densities of w 190 
^Jy and « 240 //Jy, respectively. In the initial observation of 
2MASS 0036+18 the persistent flux density was only w 130 
^Jy. The detection of persistent emission from BRI 0021 is 
marginal, with an average value of 40± 13 /iJy at 8.46 GHz 
from the combined observations is January and May. At 4.86 
GHz the flux from this source is 40 ± 18 /iJy, and it is not clear 
whether persistent emission is detected. 

In addition to the change in persistent flux between the two 
observations of 

2MASS 0036+18, the persistent emission from this source and 
TVLM 513 appears to vary within each observation. In the case 
of 2MASS 0036+18 the source brightened over the last 60 min 
of observation on Oct 9, possibly indicating a second strong 
flare (Figure ||). The sharp decrease in flux over the first 30 min 
of the observation could be the exponential tail of an earlier 
strong flare. An increase in flux was also observed towards the 
end of the observation from Sep 23 (Figure In both cases, 
the rise time is similar to the observed flare from Oct 9, indi- 
cating that perhaps strong flares from this source are common, 
and have similar profiles. 

It is therefore possible that the persistent emission is simply 
a superposition of several, possibly weaker, flares. If this is 
the case, the flare duty cycles for these two sources approach 
100%. This is similar to LHS 2065, in which Martin & Ardila 
(2001) find that weak Ha flares may be more common than 
strong flares (see §^. The implications of such sustained flar- 
ing emission are discussed in §^. 

We find no evidence for flares from the other nine sources. 
However, for most of the sources this is not surprising given 
the observed duty cycles of strong flares, and the required lu- 
minosities for a significant detection. The most luminous flare 
in the sample, from TVLM 5 13, can only be detected out to a 



distance of sa 20 pc with significance ^ 5(7. Five of the nine 
undetected sources lie at or beyond this distance. For the other 
four sources the flares could be weaker, undetected in the short 
observations (« 2 hr), or simply do not exist. Therefore, the 
lack of radio emission does not significantly constrain the pro- 
duction of flares in these sources, and the upper limits only pro- 
vide constraints on the persistent emission. 

The strong quiescent emission in TVLM 5 13 and 
2MASS 0036+18 can only be detected out to a distance of ^ 15 
pc in similar observations. Thus, some of the non-detections do 
not provide significant constraints given the detected luminosi- 
ties. In Figures || and^, and §^ we show that most of the upper 
limits on the ratio of radio luminosity to bolometric luminosity 
are consistent with the persistent luminosities of the detected 
sources. The most notable exception is the upper limit on emis- 
sion from LHS 2065. We see from Figures]^ and ^ that the 
upper limit is comparable to the detection of quiescent emis- 
sion from LP944-20 and the possible quiescent emission from 
BRI 0021. This non-detection also contrasts with the detection 
of Ha emission from this source (Martin & Ardila 2001). In 
addition, a comparison of the upper limits from Krishnamurthi, 
Leto & Linsky (1999) to the detections in this survey provides 
significant constraints on the lack of emission from their M7- 
M9 objects. As shown in Figures || and ^, their upper limits 
are up to an order of magnitude fainter in Liad/iboi than the 
detected sources in this survey. We discuss the implications of 
these non-detections in detail in §|[ 

5. ANALYSIS OF THE RADIO EMISSION 

To derive the properties of the flares and persistent emission, 
and the physical conditions that give rise to this radio emission, 
we model the Hghtcurves of TVLM 5 13 and BRI 0021 with a 
Gaussian profile: 



^v(f) = -?V,q + ^'i.,oexp 



1 ft 

2 \ at 



(1) 



Here, /v.q is the quiescent component, fo is the flare peak time 
relative to the start of the observation, ai is the HWHM du- 
ration of the flare, iv o is the peak flux, and we absorbed the 
normalization factor of l/y^vrof into the definition of Fi^_o- 
The Hghtcurve of 2MASS 0036+18 from Oct 9 is not weil- 
described by this model, and we instead use an exponential pro- 
file, in which the Gaussian in Equation [l| is replaced by rising 
and decaying exponentials with timescales and Td, respec- 
tively (Figure [3]). In Table || we summarize the derived param- 
eters for each source. We note that the profile of the flare from 
2MASS 0036+18 is somewhat unusual since flares usually ex- 
hibit a rapid rise and exponential decay as opposed to the slow 
rise and fast decay in this case (Figure H). A similar profile was 
found for two of the three flares from LP 944-20 (BOl). 

The fits to the Hghtcurves of TVLM 5 13 (Figure ^ and 
2MASS 0036+18 (Figure^, have unsatisfactory values, ~ 3 
per degree of freedom. This is due to significant variability in 
the persistent component, indicating that it is possibly a super- 
position of multiple flares, rather than a truly constant compo- 
nent. 

The statistical significance of the flare detections are deter- 
mined by the level at which we can rule out the hypothesis 
that /v = in Equation |l], coupled with the number of time 
bins in each Hghtcurve. For the flares from TVLM 5 13 and 
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2MASS 0036+18 the detections are secure, with significance 
levels of ~ 2Qa. The flare from BRl 0021 has a lower statistical 
significance of 4.5(7, but this is still high enough for a signifi- 
cant detection. We further tested the significance of this flare by 
constructing equivalent lightcurves for random positions on the 
map. We find an average fluctuation of 2a, and no fluctuations 
larger than 2.5<t, indicating that the increase in flux in BRl 0021 
is a likely flare. 

It is important to note that even if the short-timescale in- 
crease in flux in the BRl 0021 lightcurve is simply a 4.5(T sta- 
tistical fluctuation in the lightcurve, the flux density of this 
source from the combined observations on Jan 25 and May 21 
is 83 ± 15 /iJy, and hence significant at the 5.5a level. There- 
fore, whether this detection is of persistent emission or a faint 
flare, this source is radio active. 

The peak-flare luminosities are in the range Li,,f w (6- 13) x 
10'^ erg sec"' Hz"', more luminous than the brightest Solar 
flares, L^^ f w 2 x 10'^ erg sec"' Hz"' (Bastian, Benz & Gary 
1998), but somewhat weaker than flares from active early M 
dwarfs (e.g. White, Jackson & Kundu 1989). The ratio of 
flaring radio to bolometric luminosity is higher by a factor of 
~ 10^ than for Solar flares, but is very similar to flares on early 
M dwarfs (see Figures || and ||) since the lower luminosity of 
the flares relative to those from early M dwarfs is compensated 
by the lower bolometric luminosity of the late M and L dwarfs 
in this sample. The persistent luminosities are lower by factors 
of 3, 5, and 9 than the flare luminosities in 2MASS 0036+18, 
TVLM513, and BRl 0021, respectively. 

The radio band energy release in the flares is Er ^ I A x 
10^^d^^VGHzF,y.mSyCr,Mn ~ (2-10) X lO^*" erg, with the lowest 
energy release in the flare from BRl 0021, and the highest in 
the flare from TVLM513. In making this estimate we have as- 
sumed that the flare spectral energy distribution peaks near 8.5 
GHz, and has steep spectral slopes below and above this fre- 
quency, so that the bulk of the energy release is at ~ 8.5 GHz. 
This assumption is justified since for synchrotron emission (see 
§ |5.2| ) the spectral slope below the peak (self-absorbed portion) 
is /v oc i'^'^, while above the peak typical values ranges from 
- i/"'^ to i^"'* (e.g. Dulk & Marsh 1982). The high fraction of 
circular polarization in the two bright flares from TVLM513 
and 2MASP 0345+25 indicates that the observing frequency 
(8.46 GHz) lies at or above the peak frequency, while obser- 
vations of flares from several early M dwarfs (Leto et al. 2000) 
show that in these sources the flux drops significantly above 8.5 
GHz. In fact, for the M5.5 dwarf UV Get the spectral slope 
above 8.5 GHz is constrained to be lower than -2.3. Therefore, 
it is likely that most of the energy in the observed flares in this 
survey is released at 8.5 GHz. 

5.1. Violation of the Guedel-Benz relations 

One of the puzzling facts about the radio emission from 
LP 944-20 (BOl) is that despite its spectral type (M9) it vio- 
lates the Guedel-Benz relations (which hold for spectral types 
<M7) by four orders of magnitude. Here we examine whether 
any of the sources detected in this survey show a similar viola- 
tion of these relations. 

Based on these relations, and the flare duty cycle we find 
that the expected X-ray luminosity of BRl 0021 is Lx.pied ~ 
2.6 X 10^^ erg sec"', dominated by the persistent component. 
This value is ^ 10^^ times larger than an upper limit of Lx < 



2.6 X 10^^ erg sec"' from a 63.2 ksec ROSAT observation 
(Neuhauser et al. 1999). Even if we neglect the persistent emis- 
sion, the predicted luminosity due to flares is a factor of > 150 
too high. Finally, if the observed increase in the radio flux is 
a statistical fluctuation rather than a genuine flare, and the flux 
of the source is steady at 83 ± 15 //Jy, this emission violates 
the Guedel-Benz relations by a factor of > 1700. Therefore, 
regardless of whether the radio emission from this source is 
flaring or persistent (or both), it violates the Guedel-Benz re- 
lations by several orders of magnitude. This is puzzling since, 
as in the case of LP 944-20, BRl 0021 has spectral type M9.5, 
indicating that radio and X-ray emission become uncorrected 
over a narrow range of spectral type. 

As far as we know, TVLM513 and 2MASS 0036+18 have 
not been observed in X-rays, so a similar analysis is not possi- 
ble. However, we can still show that the radio emission from 
these sources is much stronger relative to the bolometric lumi- 
nosity than in early-mid M dwarfs, based on a comparison of 
the predicted peak X-ray luminosity and the bolometric lumi- 
nosity. In the case of TVLM 513, the predicted peak X-ray lu- 
minosity is w 4 X 10^^ erg sec"' , which is approximately 25% of 
the bolometric luminosity of this source. In 2MASS 0036+18 
this ratio is even higher, ^ 50%. Liebert et al. (1999), among 
others, have found a similar result for Ha flares, and noted that 
these flares may release more energy relative to the bolometric 
luminosity than in earlier spectral types. 

Therefore, while the Guedel-Benz relations hold for spec- 
tral types <M7, there is some indication from LP 944-20 and 
BRl 0021 that these relations break down around spectral type 
M9. This possibly indicates a sudden change in the mechanism 
that gives rise to the radio emission. We note that radio and X- 
ray observations of a larger sample are required to test this idea 
comprehensively. 

5.2. Magnetic Fields and Coronal Densities 

So far, we have analyzed the radio emission in a relatively 
model-independent way. However, the next step is to determine 
the emission mechanism and physical conditions that gave rise 
to the flares and persistent emission. The first step in this anal- 
ysis is to calculate the brightness temperatures: 

Tb « 2 X lOX.^jyt/pVcHzC'R/'Rj)"'- (2) 

Assuming that the emission is coronal, we can use R ^ {2- 
4)Rs ^ (2-4)R] as the size of the corona in units of Jupiter 
radii (Linsky & Gary 1983; Dorman, Nelson, & Chau 1989; 
Burrows, Hubbard, & Lunine 1989; Chabrier et al. 2000; Leto 
et al. 2000). These size estimates lead to brightness temper- 
atures of w 10** - 10^ K during the flares, and w 10^ - 10** K 
for the persistent emission. Alternatively, it is possible that the 
flares arise from much smaller regions (e.g. coronal loops), in 
which case the brightness temperatures would exceed 10" K, if 
^loop ^0-lRs ^O.lRj. While we have no data to constrain the 
size of the coronae in these sources, it is likely that the actual 
brightness temperatures are within the range defined by these 
estimates. 

The high brightness temperatures indicate that the emission 
is non-thermal. Most likely it arises from the incoherent gy- 
rosynchrotron or synchrotron processes (DuUc & Marsh 1982; 
Dulk 1985), but coherent emission is also applicable if the flares 
actually have Ti, lO" K. The high degree of circular polariza- 
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tion in the flares from TVLM513 and 2MASS 0036+18 possi- 
bly favors coherent emission, but can also be achieved in gy- 
rosynchrotron emission (Dulk 1987). On the other hand, the 
high frequency of the observed emission, = 8.46 GHz, favors 
incoherent emission since coherent emission is limited to <J 10 
GHz (Bastian, Benz & Gary 1998; Stepanov et al. 2001). Since 
it is not clear which emission process produces the observed 
radio emission, we estimate the magnetic field strengths and 
coronal densities using both. 

To accurately calculate the magnetic field strengths and coro- 
nal densities in the radio active sources using the gyrosyn- 
chrotron formulation, requires detailed spectral and geometrical 
information, since the results depend sensitively on the high- 
frequency spectral slope, the viewing angle, and the geometry 
of the emission region (Dulk & Marsh 1982; Dulk 1985). For 
approximate values we can instead use the synchrotron formu- 
lation, along with the approximation /circ ~ 3/7,Tiin, where 7min 
is the minimum Lorentz factor of the radio-emitting electron 
distribution. For TVLM513 and 2MASS 0036+18, the Lorentz 
factor is 7min « 5, while for BR10021 it is 7min ~ 10. The re- 
sulting magnetic field strengths are Z? « 51 v^^ GHz7min ~ 20 G 
in TVLM513 and 2MASS 0036+18, and w 5 G in BR10021. 
These values should be treated with caution, and they can vary 
upward by an order of magnitude or more. In fact, if we choose 
characteristic values for the radio spectral index and the view- 
ing angle, we find gyrosynchrotron values of Z? ^ 350 G for 
TVLM 5 1 3 and 2MASS 0036+1 8, and B - 50 G for BRl 002 1 . 

For comparison, the Solar magnetic field as inferred from a 
multi-wavelength study of three Solar flares (Kundu et al. 200 1 ) 
ranges from 300-900 G at the loop footpoints, while fields 
in active early M dwarfs reach strengths of a few kG (Saar & 
Linsky 1985; Haisch, Strong & Rodono 1991; Johns-KruU & 
Valenti 1996; Stepanov et al. 2001). 

The total number of electrons in the synchrotron formula- 
tion is easily calculated using « 10^^^^^^/Va,^lJy7mi^mll ~ 
(5-20) X lO^-'. However, to calculate the densities we have 
to know the geometry of the flare region. In the gyrosyn- 
chrotron formulation, using the magnetic field values quoted 
above, we find column densities of ^ 10'^ cm"^ for TVLM513 
and 2MASS0036+18, and - 10^^ cm'^ for BRl 0021. These 
densities are similar to those calculated for Solar and M dwarf 
flares (Bruggmann & Magun 1990; Stepanov et al. 2001). 

Alternatively, if the flare brightness temperatures are of the 
order of 10" K, then we must use coherent emission to estimate 
the magnetic field strengths and the densities. There are two 
main coherent processes, electron cyclotron maser (ECM), and 
plasma emission, which operate at the fundamental frequencies 
Vc ~ 2.8B MHz and i^p « 9n^^^ kHz, respectively. Recently, 
Stepanov et al. (2001) have shown for the flare star AD Leonis 
that emission due to the ECM process is strongly damped by 
thermal electrons in the corona. They conclude that for typical 
coronal temperatures and densities plasma emission is the more 
likely process. 

Plasma radiation is primarily emitted at v k, Vp. Thus, we 
can simply use Vp ~ 8.5 GHz to find n ^ 10^^ cm"^, and ap- 
proximate Vc <j O.Svp to find B <^ 1500 G. These values are 
consistent with the gyrosynchrotron-derived values. Therefore, 
we conclude that regardless of the exact emission mechanism 
or magnetic field configuration, the observations are consistent 
with magnetic fields in the range 10-10^ G, and coronal den- 



sities of the order of lO'^ cm ^. We discuss the implications of 
the inferred magnetic field values and coronal densities in more 
detail in §|. 

5.3. Heating Mechanisms 

We next investigate whether our radio data support a partic- 
ular model of coronal particle acceleration. Heating can be due 
to acoustic waves generated by the same turbulent convection 
that presumably gives rise to the magnetic fields (Bohn 1984; 
Ulmschneider, Theurer, & Musielak 1996). Ulmschneider et 
al. (1996) provide an estimate of the total flux due to acous- 
tic heating, Ff^, as a function of effective temperature and sur- 
face gravity based on an extended Kolmogorov turbulent energy 
spectrum. We use their Table 1 to estimate log(LA/i'boi) '--^ -13 
for our survey sources. In comparison, the values for the ob- 
served radio flares are log(Li-ad,f/iboi) ^ -6 (Figures |and||). 
Therefore acoustic heating cannot supply the necessary energy 
to power the flares, and similarly, cannot explain the persistent 
emission with log(L,ad,q/iboi) ~ -7. 

A more plausible heating mechanism is magnetic reconnec- 
tion (Bastian, Benz & Gary 1998; Sturrock 1999). This model 
can also explain the persistent emission if it is a superposition 
of weaker flares. In the context of this model, the flares are pro- 
duced when coronal magnetic loops reconnect, release energy, 
and create a current sheet along which ambient electrons are 
accelerated. Recent work on three-dimensional reconnection 
shows that the particles are accelerated into a power-law en- 
ergy distribution (Schopper, Birk & Lesch 1999), as required 
for the synchrotron and gyrosynchrotron discussion in §5.2. 
The accelerated electrons drive an outflow of hot plasma into 
the corona as they interact with, and heat chromospheric ma- 
terial (i.e. chromospheric evaporation). The interaction of the 
outflowing plasma with the same electons produces X-ray emis- 
sion via the bremsstrahlung process (Neupert 1968; Hawley et 
al. 1995; Guedel et al. 1996). This so-caUed Neupert effect 
indicates that there is a causal connection between particle ac- 
celeration, which is the source of radio emission, and plasma 
heating, which results in X-ray emission. This connection ex- 
plains the phenomenological Guedel-Benz relations. 

It is possible that the violation of the Guedel-Benz relations 
in the emission from LP 944-20 and BRl 0021 indicates that the 
Neupert effect no longer holds in late M dwarfs. A direct test 
of this possibility requires simultaneous observations of such 
sources in the radio and X-rays. Future surveys of activity in 
late M and L dwarfs would be considerably more effective if 
they followed such a strategy. 

It is also possible that the energy source for the radio emis- 
sion in the active survey sources and LP 944-20 is an altogether 
different process, in which case a causal relation between radio 
emission and activity in other bands would not be required. At 
present it is not possible to distinguish between these different 
possibilities, but future observations with wide spectral cover- 
age will likely provide an answer 

6. DISCUSSION 

Before we proceed to interpret the results of this survey in 
the context of correlations between radio activity and physi- 
cal source parameters, we summarize the main observational 
results: (i) the observed flares have approximately equal lumi- 
nosity, ratio of radio to bolometric luminosity, and released en- 
ergy, (ii) the persistent emission has approximately the same 
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luminosity, and it appears to be variable, (iii) LHS 2065 has 
an upper limit on the ratio of radio to bolometric luminosity 
which is lower than the detected persistent emission, despite 
being active in Ha, (iv) the emission from BRI0021, as well 
as LP 944-20, violates the Guedel-Benz relations by several or- 
ders of magnitude, (v) the inferred magnetic field strengths are 
~ 10- 10^ G, and the coronal densities are ~ lO'^ cm~^, and 
(vi) the heating mechanism is probably magnetic reconnection, 
but it is possible that the Neupert effect is violated. 

A comparison of these results to the radio properties of 
LP 944-20 (BOl) is illustrative. The peak luminosity of the 
flares from all four sources varies by only a factor of two, and 
the energy release only varies by a factor of five. The ratio of 
the flare peak luminosity to the bolometric luminosity of each 
source is log(Liad f/Lboi) ~ -6, and the ratio for the quiescent 
luminosities is log(Lrad q/i^boi) ~ -7.5 to -6.5 (Figures ^ and 

!)■ 

We find similar results when we compare this sample to other 
M dwarfs that have been observed in the radio. For exam- 
ple, flares have been detected from the dM3.5e dwarf Ad Leo 
(Stepanov et al. 2001), the dM4e dwarf Rst 137B (Lim 1993), 
and the dM5.5e dwarf Proxima Centauri (Lim, White & Slee 
1996), with ratios of Z^^ad/Lboi that are in the range of ~ -6 to 
-8. These values are plotted in Figures I and |i, and we can 
see that with the exception of the flare from Proxima Centauri, 
which is fainter by approximately two orders of magnitude, 
these values are similar to those observed in this survey and in 
LP 944-20. The flat distribution of Liad/iboi for the flares from 
the sources with spectral types M8-L3.5, and perhaps M3.5- 
L3.5 (with the exception of Proxima Centauri), possibly points 
to a "saturation" effect, in which Lmd.f/i'boi reaches a maximum 
value that is independent of spectral type or other physical prop- 
erties. It is beyond the scope of this paper and data to assess this 
possibility in detail. 

The pattern of persistent radio emission in early M dwarfs 
through L dwarfs is more puzzling. In Figures ^ and ^ we plot 
the ratios of persistent radio to bolometric luminosities for the 
detections and upper limits in this survey, the upper limits from 
Krishnamurthi, Leto & Linsky (1999), and detections in the 
range M3-M6.5 (Cash, Chai'les & Johnson 1980; Pallavicini, 
Wilosn & Lang 1985; Leto et al. 2000). Despite the small 
sample, there is a clear increase in the ratio of persistent radio 
to bolometric luminosity over the range M3-L3.5 (dashed and 
dotted lines in Figure^). This increase is in direct contradic- 
tion with the steep decline in Ha activity beyond spectral type 
M7 (Gizis et al. 2000). In fact, while the Ha activity drops by 
approximately 1.5 orders of magnitude between M7 and L3.5, 
the radio activity appears to increase by a similar amount in 
the same spectral range. This indicates that the radio and Ha 
activity are probably not correlated at the bottom of the main 
sequence. 

Another interesting result from Figure ^ is that there are 
several upper limits, which appear to fall below the appar- 
ent steady increase in Lrad/i'boi- These upper limits are for 
2MASS 1507-16 (L5) and SDSS 1346 (T6) from this survey, 
and five sources from Krishnamurthi, Leto & Linsky (1999): 
VB8 (M7), VBIO (M8), GL569B (M8.5), LHS 2065 (M9), 
and Gl 229B (T); we use the upper limit for LHS 2065 from Kr- 
ishnamurthi, Leto & Linsky (1999) since it is approximately 3 
times lower than the value measured in our survey. While some 
variance in activity levels within the sample of M and L dwarfs 



is not unexpected, it is interesting to investigate whether these 
source have some property which differentiates them from the 
detected sources. In Figure ^ we plot the ratio of persistent 
radio to bolometric luminosity as a function of the source rota- 
tional velocity, and we find that of these seven sources the three 
with measured rotational velocities have vsin/ < 10 km sec~\ 
while the detected sources have vsin/ > 30 km sec~^ . In fact, 
the source with the most constraining upper limit (VB 8) has the 
lowest rotational velocity in the sample, while the source with 
the highest value of Li-ad/iboi (TVLM513) has the highest ro- 
tational velocity. The upper limit on the radio activity from the 
fast rotator Kelu 1 (Krishnamurthi, Leto & Linsky 1999), with 
vsin/ w 60 km sec"' (Basri et al. 2000) is three times higher than 
the detected fast rotators, and therefore does not constrain the 
possible relation between radio activity and rotational velocity. 

We note that while these sources provide tantalizing evidence 
for a relation between rapid rotation and radio activity, the sam- 
ple is still very small. Future observations of a larger sample 
will provide a better test of this possible relation. However, it 
is already possible to rule out the idea that high rotational ve- 
locities quench the magnetic dynamo through supersaturation 
(Randich 1998), as was inferred from the drop in persistent Ha 
and X-ray activity in late M and L dwarfs. 

The most puzzling implication of the observed radio emis- 
sion, regardless of its possible correlation with fast rotation, is 
that late M and L dwarfs posses appreciable and sustained mag- 
netic fields, as well as high particle densities in at least some 
regions of their coronae. This result seems to contradict re- 
cent calculations by Mohanty et al. (2002) which indicate that 
the generation and propagation of magnetic energy decreases 
with decreasing effective temperature, resulting in a decrease 
of coronal activity in late spectral types. 

In addition, the presence of magnetic fields and ionized ma- 
terial indicates that magnetic braking should have slowed these 
sources down considerably. Thus, it is not clear how these 
sources maintain such high rotational velocities. 

No other correlations are apparent from the data. For exam- 
ple, the upper limit on the young object CRBR 15 is approx- 
imately five times higher than the faintest detected persistent 
emission, and therefore precludes any conclusion on the role 
of age in the production of radio emission. Similarly, the up- 
per limits on emission from the confirmed binary systems DE- 
NIS 1228 and GD 165B are approximately forty times higher 
than the detected sources, and it is therefore not possible to as- 
sess the role of binary interaction in the production of radio 
emission. 

7. CONCLUSIONS 

We detected flaring and persistent emission at 8.46 GHz 
from the M and L dwarfs TVLM513, 2MASS 0036H-18, and 
BRI0021, but no emission from the nine other survey sources. 
These observations lead to several interesting conclusions, 
which are puzzling given the decrease in persistent Ha and X- 
ray activity in late M and L dwarfs. 

The flares and persistent emission from the three survey 
sources, as well as LP 944-20 (BOl) have similar properties. In 
particular, the emission from these sources has similar Li-ad/iboi 
values, and energy release. This probably points to a common 
emission mechanism and physical conditions (e.g. magnetic 
fields) in these sources. 

In addition, while the dependence of Ha activity on rota- 
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tion appears to break down in early-mid M dwarfs, the radio 
active sources for which the rotational velocity is known are 
rapid rotators, with the most luminous source in the sample, 
TVLM513, having the highest rotational velocity, vsin/ w 60 
km sec"' . At the same time, the limits on radio activity in three 
slow rotators from Krishnamurthi, Leto & Linsky (1999) indi- 
cate that they are less active than the detected sources in this 
sample, pointing to a possible relation between rotation and ra- 
dio activity. Regardless of whether a correlation actually exists, 
the radio emission requires magnetic fields and coronal densi- 
ties that should have resulted in significant slow-down due to 
magnetic braking. It is still not clear how these sources main- 
tain their rapid rotation. 

Possibly related to the differences in radio activity and other 
activity indicators is the violation of the Guedel-Benz relations 
in BRI 0021 and LP 944-20. In these two sources, the predicted 
X-ray emission is over-luminous relative to X-ray observations. 
This indicates that radio and X-ray activity are no longer corre- 
lated in some late M dwarfs, in the same way that radio and Ha 
activity appear to be uncorrected. Moreover, it appears that the 
breakdown in the Guedel-Benz relations occurs over a narrow 
spectral type range, M7-M9. 

The observed persistent emission is clearly variable, possi- 
bly on a timescale of the order of one hour (Figures |l]and|). It 
is possible that this emission is simply a superposition of weak 
flares, in which case the activity duty cycle is ~ 100%. How- 



ever, if this component is truly persistent over long timescales, 
it means that our understanding of magnetic field generation 
in very low mass stars is far from complete. Longer, multi- 
frequency observations are required for a thorough understand- 
ing of the origin of this emission component. 

To address the possible correlation between rotation and ra- 
dio activity, as well as the differences between the activity indi- 
cators (i.e. radio. Ha, and X-rays), and to gain a deeper under- 
standing of physical processes in dwarf stars it is necessary to 
carry out a comprehensive, multi-wavelength survey of nearby 
M and L dwarfs. In particular, it is crucial to observe a wide 
variety of objects for each spectral type, and to reach ratios of 
log(Lrad/i'boi) < "8. Dcspitc the complexity of such an effort, 
the results will impact a wide range of fields from dynamo the- 
ory to the formation of brown dwarfs. 
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Table 1 

Properties of the Survey Sources 



Source 


Spect. Type. 


Dist. 


Tcff 


Age 


I'sin; 


LogCLboi/i-©) 


Log(LHa/ibol) 


Log(Lx) 


Pred. F^,R 


Flares? 






(pc) 


(K) 


(Myr) 


(km/sec) 






(erg/sec) 


(M-fy) 




CRBR 15 


M5 


160 


2930 


<3 




-1.3 




< 28.0 


<0.1 


no 


LHS 2243 


M8 


17 


2880 




<5 




-3.6 






yes 


TVLM513 


M8.5 


10.5 


2500/2880 




60 


-3.6 








no 


LHS 2065 


M9 


8.5 


2100/2630 


> 500 


9 


-3.5 


-4 


26.6 


1.5 


yes 


BRI0021 


M9.5 


12 


2300 




42 ±8 


-3.5 


-6.3 


<25.4 


<0.05 


yes 


PC 0025 


M9.5 


62 


2000 


< 1000 


13±3 


-3.7 


-3.4 


<27.2 


0.2 


no 


2MASP 0345+25 


LO 


27 


2000 


< 1000 




-3.6 


<-5.6 






no 


2MASS 0036+18 


L3.5 


9 


1800 


> 1000 












no 


GD 165B 


L4 


32 


1900 


< 1200 


38 


-4.0 


<-5.2 






no 


2MASS 1507-16 


L5 


8 


1700 


> 1000 












no 


DENIS 1228 


L5 


20 


< 1800 


< 1000 


11 


-4.3 


<-5.3 


< 27.3 


< 3 


no 


SDSS 1346 


T6 


11 


~900 






-5.3 


<-5.3 






no 


LP 944-20 


M9 


5 


2200 


~550 


30 


-3.8 


-5 


25.7 


0.5 


yes 



Note. — List of objects observed in this survey; LP 944-20 is included for completeness. The columns are (left to right), (1) source name, 
(2) spectral type, (3) distance, (4) effective temperature, (5) approximate age, (6) rotational velocity, (7) bolometric luminosity, (8) ratio of 
Ha luminosity to bolometric luminosity; for PC 0025+0447 the Ha emission is persistent, while for LHS 2243, LHS 2065, BRI0021, and 
LP 944-20 it is flaring, (9) X-ray luminosity; for LP 944-20 this is the peak-flare luminosity (Rutledge et al. 2000), (10) predicted radio flux 
densities based on the X-ray luminosities and the Guedel-Benz relations (Guedel & Benz 1993; Benz & Guedel 1994); these values are for 5 
GHz and should be higher at most by a factor of four for the observing frequency in this survey, u = 8.46 GHz, and (11) indicates whether Ha 
or X-ray flares have been detected from the source. A more detailed description of the sources, and the relevant references are given in 
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Table 2 

Very Large Array Observations 



Source 


Date 


Array Config. 


beam size 


^ on— source 




S±(T 




(UT) 




(arcsec) 


(ksec) 


(GHz) 


(/iJy) 


PC 0025+0447 


2001 Jan 25.85 


A 


0.34 X 0.24 


3.7 


8.46 


< 137 


BRI 002 1-02 14 


2001 Jan 25.93 


A 


0.33 X 0.26 


5.0 


8.46 


84 ±27 


2MASP J0345432+254023 


2001 Jan 26.01 


A 


0.30 X 0.26 


5.0 


8.46 


< 125 


DENIS-P J1228.2-1547 


2001 Jan 29.34 


A^B 


0.78 X 0.56 


4.9 


8.46 


< 154 


SDSS J134646.45-003150.4 


2001 Jan 30.43 


A^B 


0.65 X 0.37 


5.3 


8.46 


< 106 


GD 165B 


2001 Jan 30.51 


A^B 


0.64x0.31 


5.3 


8.46 


<97 


CRBR15 


2001 Jan 30.60 


A^B 


0.72 x 0.48 


4.5 


8.46 


< 149 


DENIS-P J1228.2-1547 


2001 May 29.20 


CnB 


1.6x0.86 


4.3 


8.46 


<87 


BRI 0021-0214 


2001 May 21.54 


CnB 


4.4 X 2.9 


5.5 


8.46 


83 ±18 


BRI 0021-0214 


2001 June 11.46 


CnB 


4.6x4.1 


6.4 


4.86 


<94 


2MASP J0345432+254023 


2001 June 11.58 


CnB 


3.0 X 1.7 


5.1 


8.46 


<88 


SDSS J134646.45-003150.4 


2001 June 20.24 


CnB 


4.7 X 1.9 


5.1 


8.46 


<68 


PC 0025+0447 


2001 June 20.41 


CnB 


4.8 X 1.8 


5.0 


8.46 


<75 


2MASSW J0036159+1821 10 


2001 Sep 23.44 


DnC 


12.1 x8.2 


7.8 


8.46 


135 ±14 


LHS 2065 


2001 Sep 23.55 


DnC 


12.5 X 7.5 


6.6 


8.46 


<81 


LHS 2243 


2001 Sep 23.64 


DnC 


8.8x7.7 


6.6 


8.46 


<76 


TVLM 5 13-46546 


2001 Sep 23.73 


DnC 


11.5 X 9.8 


6.6 


8.46 


308 ±16 


2MASSW J1507476-162738 


2001 Sep 23.81 


DnC 


12.3 X 7.9 


7.8 


8.46 


<58 


2MASSW J0036159+182110 


2001 Oct 9.13 


DnC 


11.2x7.8 


9.6 


8.46 


327 ±14 


LHS 2065 


2001 Oct 12.66 


D 


11.2x7.8 


3.6 


8.46 


<95 


LHS 2243 


2001 Nov 1.42 


D 


10.2x9.2 


11.1 


8.46 


<47 



Note. — The columns are (left to right), (1) Source name, (2) UT date of the start of each observation, (3) 
array configuration, (4) synthesized beam size, (5) total on-source observing time, (6) observing frequency, 
and (7) peak flux density at the position of each source, with the error given as the rms noise in the image; 
upper hmits are the flux at the position of the source plus 3a. 



Table 3 

Properties of the Radio Emission 



Parameter 


TVLM513 2MASS 0036+18 (Sep 23.44) 


2MASS 0036+18 (Oct 9.13) 


BRI0021 (May 21.54) 


Ave. Hux (fUy) 


308 ±16 


135 ±14 


327 ±14 


83 ±18 


Ave. /circ (%) 


46±4 


^6 ±11 


^3±3 


25 ±17 


Xr Const. 


21 


1.5 


7.2 


2.2 


Xr Flare 


2.7 




2.9 


1.2 




980 ±40 




720 ±40 


360 ±70 


erf"' (min) 








5 6+^" 
-■•"-i.o 


^I'.q (^Jy) 


190 ±15 




240±11 


25 ±15 


Peak/,,, (%) 


66±4 




-62 ±5 


30 ±20 


Energy (erg) 


1.1 X 10^7 




8.0 X 10^6 


2.2 X 10^6 



Note. — The rows are (top to bottom), (1) average flux density over the entire observation, (2) average fraction of 
circular polarization over the entire observation, (3) reduced x^ assuming that each lightciu-ve can be fit with only a 
constant term, (4) reduced x^ when we include a flare, (5) flare peak flux density, (6) flare FWHM duration, (7) quiescent 
flux density, (8) fraction of circular polarization at the peak of the flare, and (9) total energy release in the flare, assuming 
that the flares peaked at 8.46 GHz. 
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40 60 80 

t-Sep 23.73 UT (min) 

FIG. 1 .— Lightcurve of the 8.46 GHz emission from TVLM 513 on 2001, Sep 23.73 UT, with 20 5.5-min bins. We find a single flare 
and persistent emission, which appears to be variable. The solid line is a Gaussian model (Equation The inset shows the circularly 
polarized flux (diamonds) and the total flux (squares). The fraction of circular polarization near the peak of the flare is w 65% (see 
Table |). 
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FIG. 2 — Lightcurve of the 8.46 GHz emission from 2MASS 0036+18 on 2001, Sep 23.44 UT, with 15 8.5-min bins. We find 
only persistent emission, with a marginal indication for re-brightening over the last thirty minutes of the observation. This possibly 
indicates the rise of a strong flare. The sohd fine is a fit to a constant source. 
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FIG. 3.— Lightcurve of the 8.46 GHz emission from 2MASS 0036+18 on 2001, Oct 9.13 UT, with 25 6.5-min bins. We find a 
flare and persistent emission, which appears to be strongly variable. The steep decline in flux during the first twenty minutes of 
the observation, and the shallow rise during the last sixty minutes possibly signal two additional strong flares. The solid line is an 
exponential model. The inset shows the circularly polarized flux (diamonds) and the total flux (squares). The fraction of circular 
polarization near the peak of the flare is w -65%. The negative values indicate left-handed circular polarization (see Table |^. 
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Fig. 4.— Lightcurve of the 8.46 GHz emission from BRI0021 on 2001, May 31.54 UT, with 17 8-min bins. We find a single flare 
with a statistical significance of 4.5o-. The solid line is a Gaussian model (Equation |l]). The fraction of circular polarization near the 
peak is 30 ± 20% (see Table |). 
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LO L2 L4 
Spectral Type 



Fig. 5. — Ratio of flare luminosities (filled symbols) and persistent luminosities (open symbols and filled/open symbols) to bolometric 
luminosities vs. spectral type. Sources from this survey are designated by circles, LP 944-20 (BOl) is designated by a diamond, the 
sources from Krishnamurthi, Leto & Linsky (1999) are designated by squares, and sources in the range M3-M6.5 (Cash, Charles & 
Johnson 1980; Pallavicini, Wilosn & Lang 1985; Leto et al. 2000; Stepanov et al. 2001) are designated by stars. The upper limits 
represent 3(7 levels. When the bolometric luminosity of the source is not known, we estimated the value using other sources with 
the same spectral type (see Table |l]). While the sample size is still small, it appears that there is no decline in radio activity between 
spectral types M3-L3.5, indicating that the radio and Ha emission are possibly decoupled. 
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Fig. 6. — (a) Ratio of radio to bolometric luminosities for detected flares from M and L dwarfs. Symbols are as in Figure ||. With 
the exception of Proxima Centauri (M5.5), the radio activity appears to cluster around log(Lrad,f/i'boi) ~ -6. This possibly indicates 
a saturation effect, (b) Ratio of radio to bolometric luminosities for persistent emission from the detected sources, and upper limits 
which are similar to or lower than these detections. The dashed line is a linear fit to the detections, while the dotted line includes 
the upper limits. In both cases there is an increase in log(Liad,q/iboi) with spectral type. The same ratio for Ha emission drops 
significantly beyond M7. These observations indicate that the radio and Ha emission are probably uncorrelated at the bottom of the 
main sequence. The upper limits that violate the relation, and have measured rotational velocities, have vsin; < 10 km sec"'. This 
result is discussed in detail in §^ and Figure ^ 
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Fig. 7. — Ratios of persistent radio luminosity to bolometric luminosity for sources with measured rotational velocities. Symbols are 
as in Figure |] The sources with detected radio emission have vsin/ > 30 km sec"' , while sources which have upper limits lower than 
the detected levels have vsin/ < 10 km sec"'. This trend is based on a small number of objects, but it provides tantalizing evidence 
for a correlation between radio activity and rotation. 



